A pure Dirac's canonical analysis for four-dimensional BF theories 
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We perforin Dirac's canonical analysis for a four-dimensional BF and for a generalized four- 
dimensional BF theory depending on a connection valued in the Lie algebra of 50(3,1). This 
analysis is developed by considering the corresponding complete set of variables that define these 
theories as dynamical, and we find out the relevant symmetries, the constraints, the extended 
r^ I Hamiltonian, the extended action, gauge transformations and the counting of physical degrees of 

rr^' freedom. The results obtained are compared with other approaches found in the literature. 
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I. INTRODUCTION 
(N 

'^ , Presently, the study of topological field theories is a topic of great interest in physics. The 

>-^ . importance to study these theories arises because they have a close relationship with general 

relativity. These theories are characterized by the absence of local physical degrees of freedom 
and by the background independence [ij. Relevant examples with close symmetries with general 
relativity are the so called BF theories, which are background independent and difFeomorphisms 



JH I covariant, and were introduced as generalizations of three dimensional Chern-Simons action or as 

a zero coupling limit of Yang-Mills theories [2, |3|. We can find several examples where the BF 
theories come to be physically relevant in alternative formulations of gravity, such as Plebanski or 
Macdowell-Mansouri formulations; the former consists in to obtain General Relativity by imposing 



a. 



The later consists in 



extra constraints on a BF theory with the gauge group SO {3, 1) or SO{A) 
breaking down the symmetry of a BF theory from 50(5) to 50(4), to obtain the Palatini action 
plus the sum of the second Chern and Euler topological invariants |;5| , and since these topological 
classes have trivial local variations that do not contribute classically to the dynamics, one obtains 
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essentially general relativity 
Other interesting case, where BF theories have a close relation with physical theories is found in 
Martellini's model [3] • This model consists in expressing Yang-Mills theory as a _Bi^-like theory, and 
the BF first-order formulation is equivalent (on shell) to the usual (second-order) formulation. In 
fact, both formulations of the theory possess the same perturbative quantum properties; specifically 
the Feynman rules, the structure of one loop divergent diagrams and renormalization have been 
studied, founding an equivalence of the uu-behavior for both approaches [7|. Furthermore, other 
kind of topological BF theories are reported in [8|, where by using a generalized differential 
calculus a geometrical relation between a generalized Chern-Simons functional and a generalized 
BF theory is obtained; this version corresponds to a pure BF term, plus the second Chern class 
and a cosmological-like term quadratic in the field B. The case of gravity viewed as a generalized 
topological field theory and its close relation with the generalized BF theory is also discussed [8| . 
In this manner, with these motivations we will perform the Hamiltonian analysis for a four- 
dimensional BF theory and for the generalized BF theory introduced in [8]. It is important to 
remark that the standard way to develop the Hamiltonian study for a BF theory is considering as 
dynamical variables only those ones that occur in the Lagrangian density with temporal derivative 
[9l. llOJI (called smaller phase space context). However this approach is only convenient to perform 
provided that the theory under study presents certain simplicity; but the price to pay for developing 
the standard approach is that we can not to know the full structure of the constraints and their 
algebra, the equations of motion and the gauge transformations. Nevertheless, the approach 
developed in this work will be quite different to the standard one; this means that in agreement 
with the background independence structure that presents the theory under study, we will develop 
the Hamiltonian framework by considering all the fields occurring in the theories as dynamical ones; 
this fact will allow us to find the complete structure of the constraints, the equations of motion, 
gauge transformations, the extended action as well as the extended Hamiltonian. We able to realize 
that developing the Hamiltonian approach on a smaller phase space context, the structure obtained 
for the constraints is not right. In fact, we observe in 'l2| that the Hamiltonian constraint for 
Palatini theory does not has the required structure to form a closed algebra with all constraints; 
this problem emerges because of by working on a smaller phase space context we lose control on 
the constraints, and to obtain the correct structure sometimes they need to be fixed by hand as 
it was done in [10| for Plebanski theory. Nevertheless there are analysis on a small er p hase space 



nail er p r 
ries [11|. 



performed without complications, as for instance in Maxwell and Yang-Mills theories [11| . For these 
reasons in this paper we develop a pure Dirac method applied to models with a close relationship 
to Palatini and Plebanski theory just as the four-dimensional BF theories, and we will see that is 
not necessary to fix by hand the constraints because the method itself provides us the required 
structure. Thus in this paper we establish the bases for forthcoming works where will be applied 
the same approach to Plebanski theory. We will discuss all these details along the paper and we 
have added an appendix to clarify these ideas. 



II. A PURE DIRAC'S ANALYSIS FOR A FOUR-DIMENSIONAL BF THEORY 



In this section, we will develop an extension of the results reported in [9[. In the following lines, 
we shall study the Hamiltonian dynamics for a four-dimensional BF theory by using a pure Dirac's 
method. With the terminology "a pure Dirac's method" we mean that in concordance with the 
background independence of the theory, we will consider in the Hamiltonian framework that all the 
fields that define our theory are dynamical ones. 
So, let us start with a four-dimensional BF theory which is described by the following action [9|, |lO| 

S[A,B]^ [ B^''Ar,j(^), (1) 

Jm 

where F^"^ = dAjj + Aj^ A Akj is the curvature of the Lorentz connection 1-form A^'^ = ^^^"'da;", 
and B^' = ^B^"^ apdx" A dx^ is a set of six 5*0(3,1) valued 2-forms. Here, ^,v — 0,1,.., 3 are 
spacetime indices, x^ are the coordinates that label the points for the four-dimensional Minkowski 
manifold M and I, J = 0, 1.., 3 are internal indices that can be raised and lowered by the internal 
Lorentzian metric rjij — ( — 1, 1, 1, 1). 
The equations of motion that arises from the variation of the action are given by 

F = 0, 
DB = 0, (2) 

in this sense, both B and A are considered as dynamical variables and we will take account this fact 
for all our developments along the paper. 

To perform the Hamiltonian framework, we will suppose that the manifold M has topology S x i?, 
where S corresponds to Cauchy's surfaces and R represents an evolution parameter. By performing 
the 3-1-1 decomposition, we can write the action as 

S[A,B] = ^ f I dtd^x {iJ^oaF/jbce"""^ + e°°'^i?/,/ab (i"c - d,A'\ ~ A'^^^Al'c + A'^^.Al^) } 

(3) 
where can be identified the following Lagrangian density 

C = iiJ^^oaF/jbce"""^ + ^e""'^i?/ja;, (i"c - d,£\ - A'^^^Al^ + A'^^^Al^) . (4) 

As was commented earlier, the cornerstone of this work is to carry out the Hamiltonian analysis 
by considering as dynamical all the set of A^^'' —^2A and B^'^ ^p — J'36 variables that define the 
action. We observe at this point, that our procedure is quite different to [9|, [lO| because in that work 
the Hamiltonian analysis was performed considering as dynamical variables only those with time 
derivative occurring explicitly in the Lagrangian (in this particular scenario only the A^'^ ^ — >18 are 
considered dynamical). It is important to observe also that strictly speaking, A^'^ a and B^'^ ap are 



our set of dynamical variables and the correct form to carry out the Hamiltonian analysis is taking 
to account that set. However, because of the action is not quadratic in the field Sp/ or A^, the 
Hamiltonian study has been performed on a smaller phase space context neglecting the variables 
Aq-^ and B^'^Qi as dynamical and being identified as Lagrange multiplier [9|, |lO|. Nevertheless, it is 
not ever easy perform the Hamiltonian analysis on a smaller phase space; example of this fact is 
present in Plebaiiski's formulation [10|, where Dirac's analysis needs a treatment with more details, 
introducing new variables and fixing the structure of the first class constraints by hand. In this 
manner, we attempt with the present paper to extend the standard approach developed for BF 
theories, by performing a full canonical analysis that will be useful in order to carry out the analysis 
for Plebanski's formulation or for models found in BF gravity [211] without fixing by hand the 
constraints as reported in jlO]. 

Hence, by identifying our set of dynamical variables, a pure Dirac's method calls for the definition 
of the momenta (n"/j,H"'^/j) canonically conjugate to [A^'' a,Bai3^'^), 

The matrix elements of the Hessian 

d^C d^C d^C 

(6) 



d{d^{AJJ))d{d^{Af,iJ))' d{d,{AJJ))d{d^{Bp,iJ))' d{d,{Bp,iJ))d{d^{B^,iJ)y 

are identically zero, the rank of the Hessian is zero, thus, we expect 60 primary constraints. 
From the definition of the momenta ([5]), we identify the following 60 primary constraints 



<iPij : H^j«0, 
11 /,/ w u. 



IJ ■ ii IJ 

n'^'ij « 0, (7) 
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where we have defined e'^"-'"^ = 77"'"^. 

By neglecting terms on the frontier, the canonical Hamiltonian for BF theory is given by 



Hr = — dx 



£''oDaIl-IJ + ^v'^'^Boa'-'Futc 



(8) 



In this manner, adding the primary constraints ([7]) to the canonical Hamiltonian, the primary 
Hamiltonian is given by 

Hp = H,+ f dx' [X^'ocbij" + \'\<Pi.f + Aoa'''0°^j + Xab^r^ij] , (9) 



where A'^'^oi A^"'q, Aqq'^'^ and Xab^'^ are Lagrange multipliers enforcing the constraints. 
The non-vanishing fundamental Poisson brackets for the theory under study are given by 

{A'\{x'>,x),U''KL{x",y)} = S\^{S'KS\-5'LS'K)S^x^y), 
{i3„/^(x°, x), H^^'xlIx", 2/)} = i {S^^aS^fi - S^^pd'^o.) {5'k5'l - 5\5' k) S\x - y). (10) 



Now, we need to identify whether the theory presents secondary constraints. For this aim, we 
compute the 60 x 60 matrix whose entries are the Poisson brackets among the primary constraints 
([7|). the nonzero brackets are given by 

{r/j(a;),/VL(2/)} = -lv''''%ViKVJL-mLVJK)S'{x,y), (11) 

this matrix has rank= 36, and 24 hnearly independent nuh- vectors. This result suggests that consis- 
tency conditions imply 24 secondary constraints. From the temporal evolution of the constraints ([7]) 
and the contraction with the 24 null vectors, it follows that the following 24 secondary constraints 
arise 

A,/ = {(t>°ij{x),Hp}^0 ^ V^,j := D.n^j « 0. 
<^%,/ = {/^,/(x),ffp}«0 ^ ^°''ij:^lv'''"Fhcij^O, (12) 

and the rank allows us fix the next values for the Lagrange multipliers 

KL 
JL'IKI - ii IL'IKJl ^( 

^abc 
abc 






riJ^{rij{x),Hp}^0 ^ v'^'^^ciJ^O. (13) 

For this theory there are not third constraints. By following with the method, we need to separate 
from the primary and secondary constraints which ones correspond to first and second class. In 
order to achive this aim, we need to calculate the Poisson brackets among primary and secondary 
constraints, which the nonzero brackets are given by 

{0^j(x), ^KL(y)} = \ (n^jLviK - TTiLTiKj + W'kjvil - n"K777ij) ^'(a; - y), 

- [AkiVlj - AkjcTIli) 5'^{x - y), 
{*/j(a;),*KL(y)} = -i^{'^LjmK-'^KjriiL + '^iLr]Kj-'^iKVLj)5'^{x-y)KQ, 
{■^iAx),-^-'kl{v)} = \{-^''LjmK~-^''K.jmL + -^''iLVKj--^''iKVLj)5\x^y).^0, (14) 

Thus, we can observe that this matrix has a rank= 36 and 48 null- vectors. In this manner, we find 
that our theory presents a set of 48 first class constraints and 36 second class constraints. By using 
the contraction of the null vectors with the constraints ([7]) and ([H]), we identify the following 48 
first class constraints 

7/j : ^v'"''Fjjbc + 2D,U'''ij, (15) 



we identify the third constraint as the Gauss constraint for BF theory, corresponding to the generator 

of 5*0(3, 1) transformations. It is important to remark, that the nuU vectors obtained from ()14p 

provide us the complete form of the first class constraints, and we have not fixed by hand their 

structure. 

The rank obtained form the matrix (fT4|) yields identifying the following 36 second class constraints 

X'^'ij : W^'ij. (16) 

On the other hand, we can observe that the 48 first class constraints given in ((T5|) are not all 
independent. The reason for that is because of in virtue of Bianchi's identity DFij — one finds 
that 

Dalh - I [n'^'lKFabj" - n'^'ajFabl''] - 0. (17) 

Thus, from the 7jj=18 first class constraints we identify that [18-6]= 12 are independent. Therefore, 
we procede to calculate the physical degrees of freedom as follows; there are 120 canonical variables, 
[48-6] =42 independent first class constraints and 36 independent second class constraints. With this 
information, we conclude that four-dimensional BF theory is devoid of local degrees of freedom, 
hence in this sense we can say that the theory is topological; although this theory has global degrees 



19| 



of freedom due to the nontrivial topology of the manifold on which is defined 

Additionally we observe that the complete structure of the constraints (TTS]) and p^ as well as the full 
structure of reducibility conditions were not reported in [9|. The reason is because in [9|, [lOJ Dirac's 
canonical method was performed on a smaller phase space context, thus the complete structure of 
the constraints and reducibility conditions were not found. Of course, in our results, by considering 
the second class constraints ([TB)) as strong equations the above results are reduced to those reported 
in [9[, thus our results extend and complete those ones. The correct identification of the constraints 
is a very important step because are used to carry out the counting of the physical degrees of freedom 
and to identify the gauge transformations if there exist first class constraints. On the other hand, 
the constraints are the guideline to make the best progress for the quantization of the theory. We 
need to remember that the quantization scheme for gauge theories as Maxwell or Yang-Mills can 
not be directly applied to theories with the symmetry of covariance under diffeomorphisms (as for 
instance BF theories) because we lose relevant physical information [l3|. 
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Now, we will calculate the algebra of the constraints; smearing the constraints with test fields 

</.2 := 7°"/./ [B] = Jdx'Ba" [no^j] , 

h := 7/,/ [C] = Jdx^Cj' [DaW'jj ~ rjabc [n^ptf^gjTy^Q - n^jpH'^Qjij 

cl>i := 7%j[D]= /"dx^Doa^ 

05 := x''ij[il]= Jdx^Hj'' 

06 := x"'/,7 [G] = Jdx^Gab" [n'^'ij] , (18) 
the nonzero brackets of the constraints are given by 

{03 [C"] ,03 [C'^'^] } = Jdx' [C'^Cj/ - C'^^C^^] 7/,7 « 0, 

{03[C"],04[Do.^^]} = |dx3[C^^Do.A'^-e^^DoaK^]7%,7«0, 

{03 [C^'^] ,05 [H,^^] } = Jdx' [C^^H,,,-' - C-^'H^k'] x'ij « 0, 
{05 [H,"] , 06 [G;,^'^] } = -^rj'^''^ J dx' [n^KHGbc''"] , (19) 

where we are able to appreciate that the constraints form a set of first and second class constraints 
as expected. From the constraint algebra (fTOl) . we are able to identify the Dirac brackets for the 
theory, by observing that the matrix whose elements are only the Poisson brackets among second 
class constraints is given by 



CaB — 



1 ^-.abc 



rf'"'iViKVJL - 'r]iLr]JK)S^{x - y) 







V 



(20) 



/ 



In this manner, the Dirac bracket among two functionals A, B is expressed by 

{A{x), B{y)}D - {A{x),B{y)}p + J dudv{A{x),C{u)}C-^{u, v){CHv),B{y)}, (21) 

where {A{x),B{y)}p is the usual Poisson bracket between the functionals A,B, C"(^) = 
(x/J^jX/j"'')? with C~}{u,v) being the inverse of (|20|) which has a trivial form. It is well known 
that Dirac's bracket (PT|) will be an essential ingredient to make progress in the quantization of the 
theory [16|,|l7|- 



Furthermore, the identification of the constraints will allow us to identify the extended action. By 
using the first class constraints ([T5|) , the second class constraints (IT51) , and the Lagrange multipliers 
(IT3l) we find that the extended action takes the form 



Se [A^", n"/,7, B,J\ n^^j, Uo", Uoa", U", Ua", vj' , vj'] = J {aJ'U'^ij + Boa"U"''lJ 



(22) 



where H is linear combination of first class constraints 

'1 



^v'^'^Fbaij - iD^n'^'ij 



(23) 

and uq^"^ jUQa^"^ ,u^'^ jUa^"' jVa^"^ ,Vab^'^ are the Lagrange multipliers enforcing the first and second 
class constraints. We can observe that by considering the second class constraints as strong equa- 
tions the Hamiltonian (|23p is reduced to the Hamiltonian found in [9| where was performed the 
Hamiltonian analysis on a smaller phase space context. In this manner, we have developed in this 
work a best and complete description at classical level. 
From the extended action we can identify the extended Hamiltonian given by 

He = H- u^'^fij - u^J'f'u - u"-iij - Ua"7°"/J- (24) 

It is well know that the equations of motion obtained by means of the extended Hamiltonian, 
in general, are mathematically different from the Euler-Lagrange equations, but the difference is 
unphysical 11 1 . 



It is important to remark, that the theory under study has an extended Hamiltonian which is 
linear combination of first class constraints reflecting the general covariance of the theory, just 
as General Relativity, thus, it is not possible to construct the Schrodinger equation because the 
action of the Hamiltonian on physical states is annihilation. In Dirac's quantization of systems with 
general covariance, the restriction of our physical state is archived by demanding that the first class 
constraints in their quantum form must be satisfied, then we can use the tools of Loop Quantum 
Gravity by finding a set of quantum states for the theory as was performed in [20] using the spin 
foam models. 
We will continue this section by computing the equations of motion obtained from the extended 



action, which are expressed by 



6n°jj : Ao'' = uo", 

sn-jj-.AJ' = 



-Da { ^Ao" + u''' 



1.1 



Boa") + Ivatc [Ao'^^n'^QLv''^ ~ Ao'^^tf^QZ.Ty^Q] 



6Boa" : n""7j 



'1 



^V^'^Fbcij - 2DbIl'''>ij 



m°^,7 : Boa" 

IJ . Tjab 



IJ 

SBab" ■■n'^-ij 



UOa 



IJ 



-y'^v^jj, 



SW'^IJ -.Bab^'' - " '"■■" '^-■"'' T^. f.. IJ TJ- IJ 



Da {Ub" - Bob") - Db {Ua" - Boa' 
Uab [n^LKAo'-'rf'-Il^LKAo"^''-'] 



Vab 



IJ 



Suo" : 7"/J = 0, 

6uoa" : 7""/J - 0, 

Su'-' : jjj = 0, 

Sua" : 7""/J - 0, 

Sva" : X"/J = 0, 

Svab" ■■ X^'iJ = 0. 



(25) 



II. Gauge generator 

By following with our analysis, we need to know the gauge transformations on the phase space of 
the theory under study. For this important step, we shall use Castellani's formalism which allows 
us to define the following gauge generator in terms of the first class constraints (fT5|) 



G = / [Doeo"j°ij + Do£oa"l°''ij + e"-iij + ej'l'^'ij] dx^ 
JT. 



(26) 



10 
thus, we find that the gauge transformations on the phase space are 

6oBab — [Da^b — DbEa \ + [£ BabP' — £' BabF \ , 

<5on°/j = 0, 

<5ono^ = 0, 
^on'^'/j - - [li'^^iFe^j - n-^'jFe^/] . (27) 

We can see that, the diffeomorphisms are not present in the previous gauge transformations; however 
it is well known that BF theory is diffeomorphism covariant. Thus, the next question that arises 
is: how can we recover the diffeomorphisms symmetry from the above gauge transformations?. The 
answer for this question can be found redefining the gauge parameters as — Eq^'^ = e^'^ — —S^^'Ap^'^ 



and ej/ =^ —S.'^B^p^''. In this manner the gauge transformations (|27l) take the following form 



K" -" Ap" + C^A^" + eF^/'\ 
B'p-'-' ^ B^J-' + C^B^J-' + C [D,B^/'' + D^Bp,'' + DpB,^"] , (28) 

which correspond to diffeomorphisms. Therefore, the latter correspond to an internal symmetry 
of the theory. It is important to remark that all this information has not been reported in the 
literature with the details that we have developed, on the contrary, usually the people prefer to 
work on a smaller phase space context. Nevertheless, the price to pay for working on a smaller 
phase space is that we could not know all the relevant information of the theory, as the complete 
form of the constraints, and the complete form of the gauge transformations. In any case, given a 
theory whose symmetries we wish study one must; first, to develop a pure Dirac's method, and then 
with all the information at hand, we will able to reproduce those results that are obtained under 
the smaller phase context. All these ideas are being already applied to Plebahski theories (see 
appendix) , however the analysis is somewhat complicated but the right structure of the constraints 
will be reported in forthcoming works |16j . 

In the later section, we shall develop the Hamiltonian framework for a theory BF-likc theory that 
emerges from a generalized Chern-Simons theory. 
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III. A PURE DIRAC'S METHOD FOR A GENERALIZED FOUR-DIMENSIONAL BF 

THEORY 

For this section, the action under consideration is given by 

S [A, B] = / \f"{A) a F ij{A) + 2kB" A F/j(A) + k^B" A B/j] . (29) 

As commented in the introduction, the above action was obtained by using a generahzed differential 
calculus, taking the generalized exterior derivative to a generalized Chern-Simons form [8]. We are 
able to observe that the first term in (f29|) corresponds to second Chern-class, the second one is a 
pure BF term and the third one is identified as a cosmological-like term J9| , where fc is a constant [8| . 
The Hamiltonian study for the action (j29p has not been reported in the literature. In this manner, 
the action ([29|) is a good example for applying a pure Dirac's method just as in above section. With 
the present analysis, we will be able to identify the full symmetries of the theory; we could think, 
however, that the action ([29]) being the coupling of topological terms, then the complete theory 
is topological as well. Nevertheless, the answer is not trivial because in the literature we can find 
examples where the couplin g of topological theories is not topological any more, since there exist 
physical degrees of freedom [15|. Therefore we need to perform the Hamiltonian analysis to know 
the symmetries of the theory. For this aim, we will proceed just as in above section developing a 
pure Hamiltonian framework. 

By taking the variation of (|29|) respect to our set of dynamical variables {A,B), the equations of 
motion are given by 

D{F + kB) ^ 0, 
k{F + kB)^0, (30) 

from Bianchi's identities DF = we can see that the second equation implies the first one. Thus, 
with that learned in earlier sections we expect for this theory reducibility conditions among the 
constraints, just as in BF theory. 
By performing the 3+1 decomposition for all the terms of the action ([29l) we obtain 



F^'^(A) A F jj{A) = \e'^^''''F'\pFij^,dx^ = V^'^F^," (^Aaij - DaAoij 

B'' A Fij{A) - ^e'^^'^'^B^^pFij^^dx^ = ^v^'^Boa^Ftdj + ^v'^'^BtJ'' (ia/j - DaA^^u 
B"aBij = ^e''^^^''B'\pBij,,,dx^ = r^-'^Boa" Bbc" ■ (31) 

In this way, the action principle takes the following form 

^ [^, B] == 11 [rj^'^FtJ'' [AalJ - DaAolj) + kr^^'^Boa" F^cIJ + kT^^'^B^c" [AalJ - DaAolj) 

+ k^7]'''"'Boa''^Bbc"]dx^dt, (32) 

and the Lagrangian density is given by 

C = 7f''F,," (^Aaij - DaAoij) + krj-'-Boa'-'Ftcij + klf'^^BbJ' (AalJ - DaAoi,i) 

+ k^T^'^^^B^J'Bt,^'. (33) 
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We can see that this theory is also background independent and has the same number of dy- 
namical variables as a pure BF theory. The momenta (n"/j, n"'^/j) canonically conjugate to 
{A"a,Baf3^'^), are given by 

6AJJJ 6B^p" ^ ' 

The matrix elements of the Hessian 

d^C d^C d^C 

, (35) 



are identically zero, the rank of the Hessian is zero, thus, we expect 60 primary constraints. 
From the definition of the momenta p4p , we identify the following 60 primary constraints 

Aj : n%j«o, 

riJ ■■ n^J - 7^^'^ {FbelJ + kBbcIj) ~ 0, 



^a'' . Tiab 



H"^j « 0. (36) 



We can appreciate that, with respect to the primary constraints for a pure BF theory, the primary 
constraints ([55]) present an extra term {Fbdj) because of the presence of the second Chern class. 
By using the definition of the momenta p4p . we find that the canonical Hamiltonian takes the form 

Hc = - I dx^ [A'-^oDall'^ij + fcSoa^'^H^j] , (37) 

the canonical Hamiltonian and the addition of primary constraints allow us to identify the primary 
Hamiltonian 

Hp = H,+ f dx^ [X"ocl>ij° + X"acl>i.f + Aoa"0""/j + Xab"r'ij] ■ (38) 

The non-vanishing fundamental Poisson brackets for the theory under study are given by 

{A'\{x'>,x),n''KL{x",y)} = S\^{S'KS\-5'LS'K)S^x~y), 

{i3„/^(x", x), H^''kl(x", y)} = i {S^'^S'^p - S^'pS'^o.) {S' kS\ ~ ^' l^' k) ^'(x - y). (39) 

Just as in the above section, we need to identify if the theory presents secondary constraints. For 
this aim, we compute the 60 x 60 matrix whose entries are the Poisson brackets among the primary 
constraints ([36]), the nonzero brackets are given by 

{rzj(a:),/^;,L(2/)} = ^,y-''^(^,^^^^-^,^^^^)j3(x,y), (40) 

this matrix has rank= 36 and 24 linearly independent null- vectors. The null vectors and consistency 
conditions imply 24 secondary constraints. From the temporal evolution of the constraints ([36]) and 
the 24 null vectors arise the following 24 secondary constraints 

0°/,/ = {0"/j(a;),i/p}«O ^ V/j := i^^n^j « 0. 
0°"/,/ = {0""/,/(a;),i/p}«O ^ ^%j:=fcH^j«0, (41) 
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and the rank fix the following 36 Lagrange multipliers 

- [v^'^FbcKJViL + v'^'^FbciKVLj] Ao^^ - v'^'^kXhcij « 0, 
r''ij^{r''ij{^),Hp}^0 ^ 77'"'^Ae,,; « 0. (42) 

From consistency of secondary constraints, does not emerge more constraints. In this way, with all 
the constraints at hand we need to identify which ones correspond to first and second class. For 
this aim, we compute the Poisson brackets between the primary and secondary constraints which 
are given in the following 84x84 matrix whose nonzero brackets are given by 

{rij{x),^'''KL{y)} = ^v'''''{viLriJK-ViKVJL)SHx,y), 

{rij{x),^KL{y)} = ^ (n'^JLVlK - n'^ILVKJ + n'^KJVlL - n'^KlVLj) s^{x - y), 
{(j)''lj{x),'^''KL{y)} = kr)"'"' iviKVJL - VKJVlL)dc + (AjLcVlK - AiLcVKj) 

- {AkiVlj - AkjcVli) S^{x-y), 

1 ., 

{*/j(x),*KL(y)} = --{'^LJViK-'^KjmL + '^iLVKj-'^iKVLj)S-^{x-y)~0, 

{*/j(x),*VL(y)} = -^(n"jLr7/K-n»jL??K,/ + n"x,7?7/L-nVmj)<53(x-2/), (43) 

After long calculations, we observe that this matrix has a rank= 36 and 48 null-vectors. In this 
manner, by using the null vectors one finds the following 48 first class constraints 



I'lJ 


rrO 


7 /./ 


TjOa 
i^^ IJ, 


IIJ 


Da^'^IJ — [n" JKBabI 


rij 


ffl^j + DbU'^'u, 



— li IK^abJ 



(44) 



where the third constraint can be identified as the Gauss constraint for this generalized BF theory, 
corresponding to the generator of 5*0(3, 1) transformations. Again, we observe that by means of 
the null vectors the form of the secondary constraints has been changed and becomes to be of first 
class; this fact is important because we do not need fix by hand the constraints to convert them in 
first class. 
From the rank we can identify the following 36 second class constraints 

X^ij : U-ij-V^'^iFb.jj + kBb.jj), 
X'^'ij : U-'ij. (45) 

So, we see that does exist a clear difference among pure BF theory and this generalized theory in 
the set of the first and second class constraints. In fact, the corresponding 7'^/j constraints and 
second class constraints x"/j arc quite different, however this result is expected due to presence of 
second Chern class and the cosmological-like terms in the action. 
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An other important point to observe is that the constraints p4|) are not all independent because of 
Bianchis identity DF — 0, that now implies 

Dal^IJ - kjIJ - k [x'^'jKBabl'' " X^'lKBabj''] " [x""' I K F abj"" ^ X"' K J F abl""] = 0, (46) 

which correspond to 6 reducibility conditions for the theory. Thus, from the jfj—18 first class 
constraints we identify that [18-6]= 12 are independent, just as for a pure BF theory. Therefore, we 
are able to procedure and calculate the physical degrees of freedom as follows; the theory presents 
120 canonical variables, [48-6] =42 independent first class constraints and 36 independent second 
class constraints. This information allows us to conclude that this generalized four-dimensional BF 
theory is devoid of physical degrees of freedom and corresponds to a topological field theory as 
well. It is important to comment that this program allowed us to know the full structure of the 
constraints, and now it is straightforward to know the results that can be obtained by considering 
as dynamical variables those occurring with time derivative in the action principle. 
Now, we will calculate the algebra of the constraints. Smearing the constraints with test fields 

</-! := 7" 



h ■= 7 /./ 



7 



yV[A]=y'dx3A"[nV], 

[B]= Jdx'Ba"[lf''ij], 

■■= IIJ [C] = j dx^Cj' [DalV'^IJ - [n^^K^.w^ - n^'^KBabj'']] , 

fij [D] = fdx^-Doa" [n^j -I- DbU'"'ij] , 
("/J [H] = Jdx^Ila" [n^j ~ r^-'- {Fbcij + kBbcij)] , 
06 := x^'iJ [G] ^ Jdx'Gab" [n'^'ij] ■ (47) 

the non-zero brackets of the constraints are given by 

{03 [c"] ,03 [c^^] } ^ Jdx' [c'^a^'' - c-'^a^'] ^j., « o, 

{03[C"],04[Do.^^]} = |dx3[C^^Do.K"-C'^^DoaK^]7%,/«0, 

{03 [c^-^] ,05 [H,^^] } = Jdx' [c'^h^k' c'^h^k'] x^ij « o, 

{05 [H,"] ,06 [Gl,^^] } = -rf^kjdx^ [n^KHGbc''"] , (48) 

where we can see that the constraints (|Tf)) correspond to a set of first and second class constraints 

respectively. 

The constraint algebra (|48| allows us to identify the Dirac bracket for the theory, and we observe 

that the matrix whose elements are only the Poisson brackets among second class constraints is 

given by 



Cap — 



k^abc 
'2 



V {viLrijK - riiKriJLWix, y) 



(49) 



V / 
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In this manner, the Dirac bracket among two functionals A, B is expressed by 

{A{x),B{y)}D ^ {A{x),B{y)}p + J dudv{A{x),C{u)}C-^{u,v){C^{v),B{y)}, (50) 

where {A{x),B{y)}p is the usual Poisson bracket between the functionals A,B, C^C"") = 
(x/J^jX/j"'') with C~}{u,v) being the inverse of (l49t which is straightforward to obtain. We will 
use in future works Dirac's bracket performing a canonical quantization scheme, since in this paper 
we are only focused on a classical description of the theories under study. 

Just as in above section, with the identification of the constraints as first and second class and by 
using the Lagrange multipliers (|42l) . the constraints (|44l) and (|45|). we find that extended action has 
the following form 

Se [Ac.", n"7J, B^J\ W^^IJ, Uo", Uoa", U", Ua", vj\ Vab"] = J {aJ-'W'ij + Boa'''n"^j 

(51) 

and here H is a, linear combination of first class constraints 

H = An" [Dali^'u ~ [Tl''\jKBabl'' - n'^'/A'Sab,/^]] - B^a" [M^j + Dbli"'' u] , (52) 

and uq^'^ ,UQa" jU^"' ,Ua" jVa^"^ ,Vab^'' are the Lagrange multipliers enforcing the constraints. 
From (I48p we similarly identify the extended Hamiltonian 

He=H- Uo'''j°jj - Uoa"l°%J - u"jjj - Ua"j°''lJ. (53) 



Again, it is remarkable to observe that the extended Hamiltonian (|53|) is a linear combination of 
first class constraints reflecting the general covariance of theory, thus, we can not construct the 
Shcrodinger equation because of the action of the Hamiltonian on physical states is annihilation 

0. 
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The equations of motion obtained from the extended action are expressed by 

Sn^jj : Ao'' = uo", 

SAa" : n^j - 2Ao^" [li'^jHVFi - TI^ihVfj] ~ "^u^" [^''jhVfi ~ ^'^ihVf.i] 

+ 2Bob^" [n"-\jHllFI - ^"''iHllFj] - 2u/" [^''^7H77F/ - n-^'^iHVFj] + 2l^''''^ D ,Vm J . 

OT^j : Aj' = -Da (^o" + u'') + k {uj' - BoJ') + Va" , 

SBoa" ■■ii°''ij = [kn'^ij + Dbn^^jj], 

SH ""ij : Boa ' = Uoa , 

SBab" ■■ n''''lJ = -klJ^'^'VcIJ + [A,j"m''HJ - Aoj"n-'Hl] + [ui"n-'HJ - Uj"U'^'hi] , 

6U''''lj:Bj'' = -DaUb''' + DbUa''^ +Vab''' -[Ao^'BabF'^ -Ao^^BabF'], 

— [u BabF' — U ' BabF ] , 

6uo" : 7°/,/ - 0, 

Suoa" ■■ f''lJ - 0, 

5u'-^ : -iij = 0, 

6vab'' ■■ X°'/J = 0. (54) 



V. Gauge generator 

By following with our analysis, we need to know the gauge transformations of the theory. For this 
important step we shall use Castellani's formalism [llj, by defining the following gauge generator 
in terms of the first class constraints (l44l) 



G= f [Doeo"j°ij + Doeoa"l'"'ij + e'''jij + ej'-f^'^ij] dx\ (55) 



'E 

thus, we find the following gauge transformations on the phase space 



SoAo" 


5oAa" 


SoBoa 


SoBab 


SoU^j 


SoU'^ij 


<5onO"7 


Soll'^'ij 



Doso'', 
-Dae'' + kea", 

DoEoa , 

[DaSb — DbSa ' \ + [£ BabF — £ BabF \ , 

0, 
0, 

- [n-'iFs'^j - ^''^7F£^/] . (56) 
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We have seen that the extended Hamiltonian for this theory was hnear combination of first class 
constraints, so diffeomorphisms are the gauge transformations, however the diffeomorphisms are not 
present in above transformations. Furthermore, it is well-known that BF and Pontryagin theory 
are diffeomorphisms covariant and we expect that ((29)) being a coupled theory of topological terms 
this important symmetry is not lost. Thus, the next question that arises is; how can we recover 
diffeomorphisms symmetry from the former gauge transformations?. The answer can be found if we 
redefine now the gauge parameters as —Eq^"^ = e'"' ~ —^^Ap^'^ and ej/ = — ^''iJ^p^'^, in this manner 
the gauge transformations ([55)1 take the form 

B'^," ^ Bp," + C^BpJ-' + e [D,Bp/'' + D^Bp," + DpB.^"] , (57) 

which correspond to diffeomorphisms. Therefore, the latter are an internal symmetry of the theory. 
It is important to remark, that this analysis has not been reported in the literature with the details 
that are displayed here. Hence, we have performed a complete local study for the action ([29)) that 
will be useful to study global symmetries by using for instance, the Atiyah-Singer theorem [19|, 
then we can use the tools developed in Loop Quantum Gravity to quantize this theory; we need to 
remember that Loop Quantum Gravity is a canonical approach where diffeomorphisms covariant 
theories are quantized without perturbative methods. 



IV. CONCLUSIONS AND PROSPECTS 

The Hamiltonian analysis for four-dimensional BF theories has been performed. By considering 
the complete set of dynamical variables that define these theories, we have obtained all the 
symmetries, the constraints, gauge transformations, the counting of degrees of freedom and the 
extended Hamiltonian. For the case of a pure BF theory, the present work has extended and 
completed the results reported in [9|, where the study was performed on a smaller phase space 
context, allowing us to know the complete structure of the constraints and the algebra associated. 
Respect to the four-dimensional generalized BF theory, and despite of there are additional terms 
in the action such as the second Chern class and the cosmological-like term quadratic in B field, we 
were able to know the principal symmetries of the theory. The analysis allowed us to conclude that 
the topological structure of the theory as well as diffeomorphisms covariance was preserved. 
With the present work, we have at hand a better classical description of the theories studied, 
thus the approach developed along the paper is an alternative way to perform a pure Hamiltonian 
framework for any theory under investigation. We can see alternative approaches in jlOl . Il2l | , where 
in the former Dirac's study for Plebaiiski's and the later for Palatini theories were performed; 
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however that study was not a pure Hamiltonian analysis as the present work; for instance, in [lO| 
alternative variables in the fields were used to carry out the analysis and the final structure of 
the constraints was fixed by hand. In this sense, we expect that our approach will be an good 
alternative way to study the symmetries of Plebanski actions (see appendix below), expecting to 
obtain a better description |16i] . 

We would finish with some remarks. Topological field theories juts as the theories studied here, are 
characterized by being devoid of local degrees of freedom. That is, the theories are susceptible only 
to global degrees of freedom associated with non-trivial topologies of the manifold on which they 
are defined and the topologies of the gauge bundle [l|, |6|, ll9| . Hence, in this paper we have analyzed 
local symmetries of the theories under study, however, we would emphasize that our results have 
been useful to analyze the moduli space of a four-dimmensional BF theory for a general base 
manifold 19]. In fact, in [l9| the studio of global symmetries for a BF theory by employing as the 
main tool the Atiyah-Singer theorem has been performed with base manifolds as S"*, K3, E{n), Sd, 
etc., from which the dimension of the moduli space has been calculated, showing that there exist 
global degrees of freedom as expected. Therefore, BF theory has been characterized both globall 
and locally providing all necessary elements to make progress in the quantization; these subjects 
will be reported in forthcoming works. 
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APPENDIX A 



In this appendix we shall develop the standard way to carry out the canonical analysis to Plebanski 



theory. Plebaiiski's theory is given [21 1 



S[B, A,^]^— f {B, A F' + ^,jB' A B^) dx^, (58) 

where F is the two form strength of a SU{2) complex one-form A = A^t*, being i' the generators of 
gauge group, B — BHi is a Lie algebra valued two-form, and $ is a zero-form. Now if we decompose 
<I> into its trace and its traceless part, namely (pij — ^^^ — ^(pSij, with (j) = ^ijS"^^ which is related 
with the cosmological constant [18], we get 



S[B, A,(j)]^— f ( Si A F* + 0„B* A B^ - ]-4>B' A bA dx^ 



(59) 
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So, by performing the 3+1 decomposition in ((59)) we obtain 

(60) 

where DaXl = daXl + e'jfe^^A^ and ■q"'"' = e°"''^ a,b,c= 1,2,3. From §U\) we identify the following 
Lagrangian density 

£ = I (^Al{r,^''^B,,,) + ^,DM'''B,,,) + ii'-^'Fl^B^.a + v''''^^lBi^,Bl - ^v"""' Bi^.B.,,,^ , (61) 

remembering that in a smaller face space context, we will consider as dynamical variables those that 
occur in the action with time derivative, so the momenta P° canonically conjugate to AJ, are given 

by 

i^r-^ = |(^""i3..c), (62) 

thus, the variation respect to the fields A^, Bq^ and cpij reads 

(63) 
(64) 

(65) 

We observe from the later equations that it is necessary to eliminate the Bq^ and (j)ij variables in 
order to identify the structure of the constraints. The Bq^ and 0^ variables are Lagrange multipliers, 
however, is not easy to identified them by using the smaller phase space approach. 
On the other hand, it is important to mention that the action (1601) has a close relation with a pure 



DaP,^ ' 


a 0, 




» 0, 


B^P"'^ - ^i^Ln? ^ 


» 0. 



BF theory and the action is not quadratic in the Boi variables [Id . [18j , however the action worked 
out in [lO| is quadratic in _Boi field and one needs involve extra variables to perform the hamiltonian 
analysis, however the canonical analysis reported in [10| becomes to be so much complicated. 



Nevertheless, in our approach the canonical analysis of the first term of the action ([59)) has been 
performed, with those results and coupling the ipijB^ A B^ and (j)B^ A Bi terms in our developments, 
we will procedure by using a pure Hamiltonian framework to perform the canonical analysis of the 
action (l59l). 
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